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A REVIEW OF NASA MINITRACK 
DATA TIME -TAGGING 

William M. Rice 


ABSTRACT 

The Mlnitrack tracking stations measure an ana- 
log of direction cosines in units called Minitrack 
counts. These counts are the scaled phase differences 
between corresponding points of the wave fronts ar- 
riving at the two antenna arrays comprising the given 
baseline. The process for generating Minitrack 
counts is examined and the time for proper time- 
tagging of the count value is identified. 


V 



SUMMARY 


The purpose of this report is to interpret and document the data time- 
tagging associated with the Minitrack interferometer. Precise orbit determi- 
nations, such as those required in the current Mission and Trajectory Analysis 
Division, Goddard Space Flight Center (MTAD-GSFC) Earth physics and geodesy 
studies, need careful interpretation of the time associated with each tracking 
system measurement. 

The intent of this report is to clearly define the value of time to be asso- 
ciated with each Minitrack measurement. This is shown to be the time that the 
counter "stop" command is generated, a time easily determined from knowledge 
of the start time and the value accumulated in the counter. The time-tagging 
procedure for the current GSFC preprocessor is described, and, it is pointed 
out in general, that for data which is not excessively noisy, it agrees with the 
description contained in this report. 

Intuitively, one might expect that time-tagging is properly achieved at the 
mid-point of the counting interval. However, this report makes it clear that 
Minitrack data is properly time-tagged at the end of the counting interval, differ- 
ing from the mid-point by as much as 5 milliseconds. 
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A REVIEW OF NASA MINITRACK 
DATA TIME-TAGGING 


INTRODUCTION: 

The Minitraek tracking stations measure an analog of direction cosines in 
units called Minitrack counts. These counts arc the scaled phase differences be- 
tween corresponding points of the wave fronts arriving at the two antenna arrays 
comprising the given baseline. Figure 1 shows one of these arrays and Figure 
2 shows the antenna tiidii arrangement. The system cimnot distinguish between 
wave fronts separated by a time delay corresponding to integral wavelengths, 
hence, it measures only the decimal jiail of the true difference. The integer 
part is supiilied by the Minitrack pre-processing computer program (Reference 3) 
which derives this number from the comparison of nearly simultaneous reading 
of the coarse, medium, and fine antenna arrays. 


MINITRACK SYSTEM OPERATION 


Definition of Minitrack Counts 


The Minitrack tracking system is capable of measuring the decimal part of 
the phase difference between tracking signals arriving at opposite ends of an 
antenna base line. The system defines one Minitrack count as 10"^ wavelengths 
of the arriving signal. 

The Minitrack count is also a measurement of time, each count being 10~® 
seconds. This can be seen with the aid of the simplified block diagram shown in 
Figure 3. 

Let us assume Lh.it a signal source is located at a very great distance from 
the receiving site, hence, we may consider that the radio frequency waves are 
collimated. The separation of the antenna arrays at the ends of the base line 
can be expressed as ' l which is the number of wavelengths at the received frequency c, ^ 
This distance can also be expressed as radians, which we shall call b. 

The signal arriving from the distant source makes an angle 0 with the base 
line, and because the signal is collimated, 0 is valid at either antenna. 

At time t the signal arriving at the nearest antenna is sin (<^ j t +(^), where </> 
is an arbitrary phase angle. (Since in theory the system does not measure 
amplitude, all signal amplitudes are considered unity.) The signal at the other 
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antenna is sin j t + ^ - b cos “(t)) where t'(t) indicates that ^ is a function of 
time. The antennas are connected to the receiving equipment by equal length 
cables. Phase shift introduced by these cables, being equal on both paths, may 
be ignored. 

The local oscillator provides two signals, one differing from the other by a 
precise frequency The basic local oscillator signal is cos 
signal is supplied to two mixers, one to heterodyne the RF signal to an inter- 
mediate frequency and the second to generate a second local oscillator signal. 

In the second mixer cos (u jt) heterodynes with sin (^jt) and the sum fre- 
quency provides sin ^0^2 0^3) t to be used to translate the second RF signal to 
an intermediate frequency. The two imermediate frequencies are 


and 


sin [(ct; j - o-j) t - b CO s6 (t) + <^ ] . 


CO S [ J - a2“^^3) t + 0J. 


These frequencies are mixed to produce sin [0.-31 - b cos 6* (t)] . 


If we compare the RF signals. 


sin (a; I t - b cos6^(t) +<^) 


and 


sin (u, jt + <^.) 


with the output of the last mixer and the £^ 3 source, 

sin [o;3 t “ b cos ^ ( t )j and sin (<.J3 t ) , 
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we note that the phase difference has been preserved. Measurement of phase 
differences can be made quite readily using the latter pair of signals because 1 X 3 
may be selected for optimization of this measurement. 

The method used to measure phase is to count the cycles of a 100 KHz clock 
between corresponding phase positions of sin (^3 t ) and sin [a ~ b cos 6 ^ ( t )] . 
The frequency of is 100 Hz so that in the corresponding period of 0.01 second 
there are 1000 cycles of the clock. This, in effect, divides the RF wavelength into 
1000 parts so that we can achieve a measurement resolution of 10 “^ v/avelengths. 

The counter is also a measure of time, each count being 10'^ seconds. 
Hence, knowing the counter start time and the count time duration, which is 
readily available from the v'alue of the count, the counter stop time can be deter- 
mined by summing these time values. 

The 100 Hz and 100 KHz signals are derived from a master timing unit 
which, in turn, is synchronized with a time standard such as WWV. The counter 
is st.arted at intervals of 0.2 seconds, coinciding with the required triggering 
le^'el of sin 3 1 ). 

An alternate way of looking at the expression 


sin (:<-3 t - bcos5 (t)) 

can be realized by transforming the term b cos ^ (t). 

It had been stated thatb = 2'^A.^ , where is the number of wavelengths 
of the received frequency separating the two antennas. The wavelength of 
frequency f ^ j = 2^^ f ^ , is c/f j . Hence the number of wavelengths can also 
be expressed as f jd c, where d is the metric length of the antenna base line. 

We c m write the equations 


b = 2 r.X^ . 




L 


f^d 


and 


b cos i9 ( t ) = 2 "^ f j 


d 

— cos 6 ? ( 


t) 


G 



but d c cos'7(t) - which is the time required for the wave front to reach the 
second antenna after arriving at the first antenna. Therefore, 


1 - 277 f, . 


b cos ^(t) =aj--(t) , 


and 


s in ['^ 3 t “ b cos c?(t)j - sin [^ 3 t - ( t ) j 


It can be seen that^c^r (t ) is the phase shift introduced by frequency prop- 
agating the extra distance necessary to arrive at the second antenna. The term 
7-(t ) is a slowly varying function of time whose value can be considered a constant 
for the interval At = '^ . 

Figure 3 also shows, in a simplified diagramatic way, the various operations 
needed to obtain a count measurement. The many needed amplifiers, shaping cir- 
cuits, filters, and other circuits in the hardware implementation have been omitted 
for clarity. 


Count Generation Analysis 

Certain assumptions are made in the following analysis: 

1. Source motion for time interval T < b/c is considered negligible. 

2. Antenna cable length differences are assumed to be zero. 

3. Circuit pass bands are assumed to be broad compared to the signal rate 
of change. 

4. Signal amplitude variations and noise are neglected. 

5. Antenna field irregularities are ignored. 
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These assumptions have the following effects on the analysis: 


1. Propagation time: The Minitrack system measures the line of sight 
angle* to a signal source. It cannot measure range nor can the prepro- 
cessor determine range. Hence, we are not concerned with the separa- 
tion between the Minitrack site and the RF source. 

The analysis does assume that, during the time required for a wave front to 
move from the nearest to the farthest antenna (an interval less than 0.5 micro- 
seconds), the change in angle to the source is too small to be significant. In 
other words, the two signals used to generate the counter stop command must 
travel unequal distances. The time difference required for this distance differ- 
ence is too small to be measured by the system. 

2. Cable length inequality: The effect of the RF line length inequality causes 
a phase shift of one signal with respect to the other prior to the genera- 
tion of the stop-counter signal. This phase shift carries over directly 

to the stop-counter signal. In actual station operation, compensation 
for this effect is achieved by careful calibration of station equipment. 

3. Circuit pass bands: In the dynamics of tracking a source, the electronic 
circuits are required to follow or track the signal while introducing a 
minimum of error. Phase integrity must be maintained throughout all 
circuits. As the source motion produces doppler variations, for ex- 
ample, it is assumed that all circuits are ideal and do not introduce 
tracking errors. 

4. Signal amplitude variations and noise: Any electronic circuit variation 
that could possibly cause a phase variation in the actual equipment has 
been ignored. Such problPins, if they do exist, are not a part of this anal- 
ysis. The effects of noise of any kind are ignored for the same reason. 

5. Antenna field irregularities: Irregularities in the antenna field could 
have the effect of causing the source to appear to be angularly displaced 
from its true direction. This effect could appear as improperly time- 
tagged data. Hence, for this analysis, the antenna fields are assumed to 
be perfect. 


*The angle that Minitrack measures is a function of the signal wavelength implicit in the pre- 
processor ond the refraction corrections employed in obtaining the station constants from the 
calibration aircraft overflights. This problem is discussed in Reference 4. 
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Referring to Figure 3, the counter start condition can be established as 


sin t j) = 0 , 

and the counter stop condition as 


sin -bcos ^(tj)] =0 , 


both valid for the positive slope of the function. 

The arguments for the sine functions are, then, 


and 


Wj t j = 2 n7T , 


oj^t 2 - b cos ^ (tj ) = 2 m77 , 


where n and m are integers, is the counter start time, and t 
stop time. Let At = counting interval time, and 

tj = tj + At . 


next. 


and 


but 


and 


ojj(tj + At)-bcos0(t2) = 2m7T , 
tt> 3 t j + ojjAT - b cos 0 (tj ) = 2 tn77 , 


OJjt j = 2 077 , 


2nT7 + <^ 3 At - b cos 0 (t 2 ) = 2 m77 , 

277(n - m) + o), At 
cos0(t2> = . 


2 is the counter 
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( 1 ) 


Because b - 2nK^ and - 2-n 


cos --(t ^ ) - 



t 


In Equation 1 the quantity (fi -ti) is an integer. Its value is closely related 
to the path length difference b cos c. However, n is evaluated at time t j , and 
m at the later time t 2 - The value t 2 - t j-At may, at times, approach one period 
of f j , which must be taken into account if the significance of n - m is to be understood. 

Let us rewrite Equation 1 so that 


n-m = \j^cosfc'(t2)“f3At , 


and 


At 

n - m - cos 6 (tj) - ' 

^ 3 


where the period of f 3 is = 1/ f ^ . We see that (n-m) is the path length 
difference at time tj reduced by a factor that is the ratio of the count accumulation 
time to the period of frequency f Generally, At /T^ will not exceed 1 , but there 
are remote circumstances which will permit this value to exceed unity. For 
example, if > ^ cos i9 (tj ) is exactly an integer, then the counter will receive 
a start command simultaneously with a stop command. Depending upon the 
characteristics of the counter, one or the other commands will be dominant, 
possibly the start command. If the source doesn’t move, the next stop 
command will be simultaneous with a start command and the counting pro- 
cess will continue until a stop command is displaced from the start command. 
Fortunately the sources for which Minitrack was designed are moving and successive 
simultaneous occurrence of start and stop commands are not possible, but it is 
within the realm of possibility for a single simultaneous occurrence to exist. If 
the next stop command should occur after the start command, then At /T 3 will be 
greater than unity. Furthermore, if ^ (t) is near tt/ 2 , (n - m)will be a negative 
number. This is demonstrated in the next section called Numerical Example. 

We can see from Equation 1 that the cosine of 0 ( t ) at time 1 2 can be 
calculated from the counting time interval At. The value of At is, quite simply, 
the accumulated count multiplied by the period of the 100 KHz signal. Hence, the 
conclusion is that the Minitrack count should be time-tagged at the occurrence of 
the counter stop command. 
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An examination of Reference 3 discloses that the preprocessor used by GSFC 
comes close to time-tagging data in this way. In fact, the difference may be less 
than the computational errors due to round-off and loss of significance. The 
GSFC preprocessor compresses each data line (frame) containing 5 data points 
to a single point. The five data points are assumed to be equally spaced. A 
parabola is fitted to these points in a least squares sense and an estimate is made 
for the smoothed midpoint value. The time associated with this value is the 
counter-stop time of the mid-point reading in that frame. For various reasons, 
one of which is for ease of resolving ambiguities, the value of the count at 0.4 
seconds after the frame start time is to be obtained. The smoothed value obtained 
from the parabolic fit is located at a time of 0.4 seconds after the frame start 
time plus the counter accumulation time. Therefore, an interpolation must be 
made to obtain a count value at 0.4 seconds after the frame start. The GSFC 
preprocessor now proceeds to estimate this value. A new assumption is made 
regarding a frame of data. It is assumed that the 5 readings in any given frame 
can be fitted with a straight line passing through the smoothed value previously 
derived. The slope to be used is the same as that for a line connecting the frame's 
first and last counts, having a time base of 0.8 seconds. The compressed frame 
point is then just the smoothed value, minus the product of the slope multiplied 
by the counter accumulation time. This completes the preprocessor's time- 
tagging of the data. 

Each line of tracking data (frame) is time-tagged in the same way for both 
east -west and north-south data. For arcs of 0.8 second duration, the difference 
between parabolic fitting and linear fitting should be quite small. However, as 
the data becomes noisy, both the linear slope and the mid-point count -accumulation 
time may become adversely affected; both are used in the final estimate of the 
smoothed count. 


NUMERICAL EXAMPLE 

Let us consider a two-dimensional example in which a radiating source 
moves in a circle using a Minltrack site as a center (Figure 4). (Although this 
is physically unrealizable it provides ease of calculation and demonstrates the 
desired results.) The plane of motion is defined by using y as the local vertical, 
and X as any direction normal to the vertical. 

The sources position is then defined as 


r 2 = y2 


+ X' 
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The tracking site is located at coordinates (0.0). Let us assume that the location 
of the source from the y axis can be expressed as 


X = r sin (10"^ t) . 


( 2 ) 


Then, at t = 0 the source is at zenith. From Figure 4, the angle 6 is shown to 
be the elevation angle, whose direction cosine is x/r. Thus, 


cos^=— 1 
r 


and from (2) 


s in (10"^ t) = 2L , 
r 


y 





«X<( 


and 


-3 


cost:?=:sin(lO“ t) 


( 3 ) 


As previously discussed, 


n - m + f , 

cos 6{t) = : 


( 1 ) 


Therefore, let 


f 3 = 100 Hz 


and 


= 50 wavel engths . 


Then, 


sin (10“^ tj ) 


n - m + 100 At , 
50 


or 


'sin (lO'^tj) = 2 (l0'^(n - m) + At) . 


( 4 ) 


£ 

next, since 

t, = t,+At. 
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Thus, 


s i n [lo-^ (tj +At)j = 2(10 (n-m) + At) . 


If the count is started at t j = 0, the source being at zenith, we can solve for 
the counter elapsed time At as follows: 

sin (10“^ At ) = 2(10"^(n - m) + At ) , 


sin (10"^At ) 10'^ At 


10“^ At = 2 X 10"^ (n - m) + 2A t , 


and 


At = 


-2x 10 (n-m) 

(2 - 10 -^ ) 


(5) 


At zenith. 


n - m = 0 . 


hence. 


At = 0. 


Perhaps the counter cannot be started and^ stopped simultaneously. Let us assume 
that At ^ 0, This implies that n-m must.have a value other than zero. Let 


n-m = -1 (see count generation analysis) . 
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Then, 


-2 < ) 
(2 - 10 '^ ) 


10"^ (1 + 5 ' lO"'* + (5 - lO"'^ - (5x 10"'»)^ t 


or 


At = 1.00050C.50125 X 10" ^ , 


which corresponds to approximately 0,5 minitrack counts greater than one full 
cycle of 100 Hz, Using actual hardware, the counter, responding to integral 
counts, will overflow and have a value of zero; but for this example, we must 
assume that the counter can measure time to a finer resolution than the 10" ^ 
seconds of the Minitrack system. To verify that the direction cosine measured 
by the system is the correct value, let us evaluate Equations 3 and 4, 

From Equation 3 we find that 

cos d - sin(10"^At ^ .oiooos = 1 • 00050025 x 10“^ . 


From Equation 4 we find that 


cos 6 = 2(10“^ (n - m) + At ) , 


or 


cos e = 2(-10"^ + 1.000500250125 X 10"^) , 


which becomes 


= 1.00050025 X 10"* . 
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If, in Equation 5, we let 


n - m = +1 . 


then It becomes a negative value; however, this will never occur because of the 
Minitrack equipment configuration. As one more example let t j = lOO; then 
t 2 = tj At. Thus, 


sin [l0“^ (t j + At )] = 2(10"^(n - m) + At ) , 


expanding, we get 


sin 10”^ tj cos 10”^ At +cos 10”^ t^ sin 10”^At = 2(10“^(n-m)+At). 

Now, because 


cos 10"^ At * 1.0 , 


and 


sin 10"^ At - 10"^ At . 


we find that 


sin lO"* + 10“^ At cos 10"^ = 2(10'^ (n - m) + At ) , 


and 


using 


At 


sin 10”^ - 2 X 10”^ (n - m) 
2 - 10"^ cos 10"* 


(2 - 10"^ cos 10" 0 


^ 0.5 
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we get 


At =10'^ [ 4 . 9916708 - (n-m)] (1 • .no04975020S26) . 


Since 


0 £ At - 10"‘ . 


then 


n - m = 4 . 


and 


At = .0099216416 . 


This would be a total of 992 Minitrack counts. 

As a check, evaluating Equations 3 and 4 with these values of ( n -m ; and At 
we find that 


cos (") = s i n [lO"^ (100.0099216416)] . 


or 


ros (•') = .099^433 . 


Also, 


cos (i) r 2(l0'^(n -n) - ’.t) - 


r 8 . 10’^ - .0198433 . 
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or 


cos (U) ^ .0998433 . 


which verifies that the Minitrack counts are applicable at time t = 1 2 . 

The time-tagging can be shown in another way by considering a plot of the 
phase angles of the start signal, sin (^3!), and the stop signal, sin U3 1 - b cos 
(t) Figure 5 shows these phase angles, modulo 2 ” versus time. (The curva- 
ture in the trace of , ^ t - b cos - (t) has been exaggerated to make the time- 
tagging problem more apparent.) If the time at generation of the start signal is 
defined as .31 = 0, then the time at generation of the stop signal may be defined 
as ' 3! - b cos (t) = 0 , both modulo 2 v, Since ^3 is derived from a precision 
source, phase angle -^3t is known at any time. If the 100 kHz counter is started 
at 1 1 and stopped at tj, it has measured time interval t2 ~ t^. which is the base 
of right triangle tj, t2* P. Side Pt2 the phase difference we seek at time tj, 
which is just 3 (t 2 - t j). 

The Minitrack count that is accumulated between t2 and tj is the scaled 
measure of this time interval in which each count is ten microseconds. The 
Minitrack count is also the scaled phase angle of <^3 at time t2 (1000 counts 
equals 2 - radians). Hence, knowing the counter start time tj and the accumu- 
lated count, the phase difference at time t2 can be determined. Thus the Mini- 
track counts are properly time-tagged at occurrence of the stop pulse. 


CONCLUSIONS 

This analysis is a practical treatment of that part of the Minitrack system 
which generates the Minitrack counts. Tne analysis shows that the time-tagging 
of the counts should coincide with the occurrence of the counter stop signal. The 
•^urrent GSFC preprocessor produces a data time-tag value that is close to being 
correct; in fact, the errors arc smaller than the current system resolution. 
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